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RELATED APPLICATIONS 

[001] This patent application claims priority to U.S. Patent Application No. 09/658,188 
filed September 8, 2000, entitled Ground Penetrating Radar Array and Timing Circuit, which is 
herein incorporated by reference in its entirety and to U.S. Provisional Patent Application No 
60/284,921 filed April 20, 2001, entitled Method and Apparatus for Identifying Buried Objects 
Using Ground Penetrating Radar, which is herein incorporated by reference in its entirety. 

BACKGROUND OF THE INVENTION 

[002] The present invention relates to ground-penetrating radar (GPR) systems and 
methods, and more particularly to ground-penetrating radar systems and methods for identifying 
buried objects. 

[003] GPR is similar to upward-looking radar used for air traffic control and meteorology, 
except that the antenna array is directed toward the ground; thus presenting an entirely different set 
of technical challenges. For example, GPR is used for geophysical applications such as mapping 
subsurface strata, locating toxic waste sites for remediation, and detecting of unexploded subsurface 
ordinance. 
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[004] A GPR system comprises at least one transmitter that transmits an electromagnet; 
impulse, continuous wave, or swept frequency signal, for example, in the frequency range of 1MHz 
to 1 0GHz. The system also comprises at least one receiver that receives a reflected waveform. The 
impulse is short with respect to the wavelength being transmitted. In practice, the impulse can be, 
for example, one or two cycles. Therefore, GPR systems are sometimes referred to as "impulse" 
radars. 
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[005] Subsurface industries such as construction, utility location, environmental 
remediation, and unexploded-ordnance detection have long sought safe, reliable, cost-effective 
methods for "seeing into the ground." The utility location market suffers greatly from inadequate 
location technologies that result in hundreds of millions of dollars in damages, delays, and lost 
revenue for utility companies and contractors every year, losses that can be reduced significantly by 
use of GPR. Three utility locating market segments, can be distinguished by their accuracy and 
price: (1) One Call; (2) private locating; and (3) subsurface utility engineering (SUE). 

[006] "One Call" is a nationwide clearinghouse that provides an alert to all public and 
private utilities of when and where construction may impact their lines. By law, contractors must 
register their site with One Call, which in turn contacts all the relevant utilities so they can mark their 
utility lines. One Call locating systems are based on electromagnetic induction technology. Utility 
companies, responding to a One Call work order, guarantee accuracy on conductive lines within 
twenty-four inches horizontally on either side, with no guarantee of depth. With One Call, utility 
line locations are simply painted temporarily on the ground, easily subject to erosion or destruction. 
This poor accuracy results in broken utility lines and revenue loss. 

[007] Construction, utility, and industrial companies often rely on "private locating." 
Private locating provides a greater degree of accuracy than is delivered by One Call. These 
companies often hire a utility locating company or a geophysics company to apply more expensive 
and time-consuming locating techniques. Private locating companies typically use electromagnetic 
induction technology, GPR, and magnetometry. Often this includes excavation, the most reliable 
and expensive conventional method for determining the exact location of utilities. 
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[008] Industrial and utility companies, however, frequently require more accurate maps of 
the subsurface than One Call or private locating can provide. For instance, extra accuracy may be 
needed while excavating near an oil pipeline because of the rise of breaking a pipe. Or, it may be 
too costly to accidently cut an interstate fiber optic cable carrying important communications. In 
such situations, excavators perform a total cost/value analysis, including consideration of risk/cost 
avoidance. Often, they are more willing to pay higher fees to ensure greater accuracy. 

[009] "SUE" can provide more accuracy than One Call or private locating. SUE is a rapidly 
growingspecialty service offered by geophysical and engineering companies. It entails planning and 
designing utility arrangements before highway or other large infrastructure construction. SUE 
engineers painstakingly map all discernible utilities at a given site using a variety of conventional 
geophysical methods. SUE uses electromagnetic induction technology, GPR, and magnetometry. It 
is generally more costly than private locating services, because it uses computer aided design to 
produce a permanent record of the location of utilities. Even this premium service often only 
identifies 80% of utilities with certainty, frequently less when unknown non-conductive utilities are 
present. Further, SUE is very expensive. Accordingly, there is, therefore, a need for systems that 
overcome the limitations of the conventional techniques. 

SUMMARY OF THE INVENTION 

[0 1 0] The present invention provides a methodology for identifying buried objects using 
ground penetrating radar. Additional features and advantages of the invention will be set forth in 
the description which follows, and in part, will be apparent from the description, or may be learned 
by practicing the invention. The objectives and other advantages of the invention will be realized 
and obtained by the method and apparatus particularly pointed out in the written description and the 
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claims hereof as well as in the appended drawings. 

[Oil] To achieve these and other advantages and in accordance with the purpose of the 
invention, as embodied and broadly described, an apparatus and method according to the invention 
identifies buried objects using GPR. An apparatus for identifying a buried object using ground 
penetrating radar (GPR) in a system containing at least one GPR sensor, comprises a data processor 
for detecting spatial correlations in data received from a GPR sensor in the apparatus and an image 
processor capable of building a data structure corresponding to an image of the buried object from 
data processed by the data processor. A method for identifying a buried object using GPR in a system 
containing a GPR sensor comprising detecting spatial correlations in data received from the GPR 
sensor in the system and building a data structure corresponding to an image of the buried object 
from the received data. 

[012] Both the foregoing general description and the following detailed description are 
exemplary and explanatory and are intended to provide further explanation of the invention as 
claimed. 

DESCRIPTION OF THE FIGURES 

[0 1 3] The accompanying drawings, which are incorporated in and constitute part of the 
specification, illustrate an embodiment of the invention and together with the general description 
given above and detailed description of the embodiment given below, serve to explain the principles 
of the invention. 

[0 1 4] Fig. 1 is an embodiment of an apparatus for identifying buried objects using ground- 
penetrating radar; 

[015] Fig. 2 is a detailed view of sensor unit 1 02; 



[0 1 6] Fig. 3 is a block diagram of three suitable antenna geometries; 

[017] Fig. 4 illustrates a direction of movement in antenna array configuration 302 with a 
single line of transmitter and receiver antenna elements; 

[018] Fig. 5 is a schematic diagram of data acquisition at three different time steps tO, tl , 
and t2 using the array configuration and scan direction shown in Fig. 4; 

[019] Fig. 6 is a schematic diagram of the antenna array elements in a configuration of 
sensor unit 1 02 for monostatic data acquisition; 

[020] Fig. 7 is a schematic diagram of multistatic data acquisition; 

[02 1] Fig. 8 is a schematic diagram of the distance between a transmitting antenna element 
and a receiving antenna element; 

[022] Fig. 9 is a graph of a received radar signal time series; 

[023] Fig. 10 is a schematic of a critically refracted ray in a GPR system; 

[024] Fig. 1 1 is a schematic of a critically refracted ray and a reflected ray in a GPR system; 

[025] Fig. 12 is a schematic diagram of bistatic measurements for computing wave speed; 

[026] Fig. 13 is a diagram of a mobile vehicle with a GPR system; 

[027] Fig. 14 is a diagram of a mobile vehicle with a GPR system; and 

[028] Fig. 15 is a diagram of a portable GPR system. 
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DETAILED DESPRTP TION OF THE INVENTION 

[029] A method and system is disclosed for identifying buried objects using ground- 
penetrating radar. Fig. 1 is an embodiment of an apparatus for identifying buried objects using 
ground-penetrating radar consistent with the present invention. The apparatus of Fig. 1 includes 
sensor unit 102, data processor 104,pre-processor 106, image processor 108, feature processor 1 10, 
detector 1 12, discriminator 1 14, and display 116. 

[030] Fig. 2 is a detailed view of sensor unit 102. An embodiment of sensor unit 102 
consistent with the present invention includes transmit antenna 202 and receive antenna 204. 
Transmit antenna 202 emits GPR Waves used to probe for buried objects. Receive antenna 204 
receives the transmitted GPR waves as they are reflected from the ground and from objects located 
beneath the ground. Sensor unit 102 controls the characteristics of the transmitted GPR signal (in 
some instances receiving support from data processor 1 04 as will be descried in greater detail below) 
and receives the reflected signal. Suitable antennas for the present invention include, for example, 
small broadband antennas with a broad beam pattern having a frequence response spanning at least 
two octaves. 

[03 1 ] Those skilled in the art will recognize that there are many suitable geometries for the 
arrangement of antenna elements in transmit antenna 202 and receive antenna 204. For example, 
three suitable antenna geometries consistent with the present invention are shown in Fig. 3 . Antenna 
302 has antenna transmitters ("T") and receivers ("R") arranged in a single line, with alternating 
transmitter and receiver elements. Antenna 304 has one row of transmitter elements and a parallel 
row of receiver elements. Antenna 306 has one row of transmitter elements and a parallel row of 
receiver elements with the receiver element row offset from the row of transmitter elements. Each 
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antenna can have the same polarization or polarization can vary among the antennas. Additional 
antenna arrangements consistent with the present invention are described in U.S. Patent Provisional 
Application No. 60/152,607 and in corresponding non provisional U.S. Patent Application No. 
09/658,188. 

[032] Although Fig. 3 depicts sensor units comprised of arrays of antenna elements, an 
embodiment of sensor unit 1 02 consistent with the present invention includes a pair of transmit and 
receive antenna elements that are moved during data acquisition to simulate data acquisition using 
an array. Throughout this application wherever exemplary embodiments are discussed using arrays 
of antenna elements, once skilled in the art will recognize that such transmit/receive antenna pairs 
may be used instead of arrays, or in combination with arrays. 

[033] To acquire data over a planar surface, for example the ground, the antennas in sensor 
unit 102 are preferably moved in a direction that is perpendicular to the line formed by the transmit 
and receive antenna arrays. For example, Fig. 4 illustrates an example of a direction of movement 
(indicated by the arrow) in antenna array configuration 302 with a single line of transmitter and 
receiver antenna elements. Fig. 5 is a schematic diagram of data acquisition at three different time 
steps tO, tl, and t2 (Figs. 5(a), 5(b), and 5(c), respectively) using the array configuration and scan 
direction shown in Fig. 4. As shown in Fig. 5(a), at time tO data are acquired over all antennas in 
the array at a fixed position xO. At time 1 1 the antenna array is moved preferably by a fixed amount 
Ax to position xl (Fig. 5(b). Once at position xl, data are acquired over all antennas in the array at 
the new position. Figs. 5(b) and 5(c) include shaded antenna array elements showing the locations 
of the antenna elements at previous time steps. At time t2 the antenna array is moved again by a 
fixed amount Ax to position x2 (Fig. 5(c)). Once at position x2, data are acquired over all antennas 
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in the array at the new position. This process continues for a number of iterations sufficient to scan 
a region of interest over a buried object. Note also that scanning can also be non-linear, e.g., Ax is 
not fixed at each time interval. 

[034] Data processor 1 04 can generate control signals that ensure proper operation of sensor 
unit 102 when acquiring data as illustrated in Figs. 4 and 5. Moreover, data processor 104 processes 
data received by sensor unit 102. Two techniques for acquiring and processing data using sensor 
unit 102 and data processor 104 consistent with the present invention are monostatic and multistatic 
processing. In monostatic data acquisition, data processor 104 assumes that a transmitting and 
receiving antenna are paired such that the data acquired simulates a single pair of antennas moved 
in unison. Fig. 6 is a schematic diagram of the antenna array elements in a configuration of sensor 
unit 102 for monostatic data acquisition. As shown in Fig. 6(a), during period tO, sensor unit 102 
fires; transmit antenna element 602a and then receive antenna element 602b records the reflected 
GPR signal. During period tl, sensor unit 102 activates the next transmit and receive antenna 
elements 604a and 604b, respectively, in the array (Fig. 6(b)). The remaining transmit and receive 
antenna elements are activated in pairs at successive time intervals, e.g., transmit element 606a and 
receive element 606b are activated during period t2 (Fig. 6(c)), transmit element 608a and receive 
element F08b are activated during period t3 (Fig. 6(d)), etc. In an alternative embodiment of a 
monostatic data acquisition consistent with the present invention, a transmitter can be paired with 
more than one receiver or vice versa. An embodiment of data processor 104 consistent with the 
present invention includes a monostatic data processor that recognizes the manner in which sensor 
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[035] Fig. 7 is a schematic diagram of multistatic data acquisition. In multistatic 
(sometimes also referred to as bistatic) data acquisition mode, data is acquired at all receivers for 
every transmitter at each time period. In an embodiment of sensor unit 102 consistent with the 
present invention, during period tO transmit element 702 fires and each of the receive elements 704- 
718 record the reflected GPR wave (Fig. 7(a)). During period tl, transmit element 720 fires and 
each of the receive elements 704-7 1 8 record the reflected GPR wave (Fig. 7(b)). During subsequent 
time periods t2 (Fig 7(c)) and t3 (Fig. 7(d)), transmit elements 722 and 724, respectively, fire and 
the resultant GPR waves are recorded by receive elements 704-718. This process continues over 
successive time intervals until data acquisition processor 102 has activated each of the transmit 
antenna elements. As stated above, note also that this technique of multistatic data acquisition can 
be performed by moving a pair of antenna elements to simulate data acquisition using an array. 

[036] An embodiment of data processor 1 04 consistent with the present invention includes 
a multistatic data processor that recognizes the manner in which sensor unit 102 acquires data 
representing the buried object. The multistatic data processor then processes data accordingly to 
extract a signal representing the buried object. An embodiment of a multistatic data acquisition 
processor consistent with the present invention has either a dedicated data acquisition channel for 
each receiver or a multiplexer for sequentially multiplexing through some lesser number of data 
acquisition channels. 

[037] Among other operations described above, including operations those skilled in the 
art will recognize are consistent with the present invention, data processor 104 in the apparatus of 
Fig. 1 performs system calibration operations including, for example, one or more of the following: 
time alignment, wave speed analysis, and impulse response function computation (which can 
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include, for example, beam pattern and spectral response computations). Each of these operation 
is described in greater detail below. Generally, data processor 104 processes data received from 
sensor unit 102 to extract a meaningful signal corresponding to buried object in a surveyed region. 
An embodiment of data processor 104 consistent with the present invention is comprised of a 
processor that executes stored computer program code designed to implement calibration operations. 
One skilled in the art will recognize that an embodiment of data processor 104 implemented entirely 
in software, entirely in hardware, or in an embodiment allocating calibration functions among 
hardware and software elements, either distributed or centralized, is consistent with the scope of the 
present invention. 

[038] An embodiment of data processor 104 consistent with the present invention can 
include a time alignment calculator that performs time alignment to compensate for the difference 
in wave propagation speed in air versus its propagation speed in geological material. 
Electromagnetic waves generated by a GPR system consistent with the present invention travel at 
the speed of light through geological objects. Generally, electromagnetic signals travel between 
controlling system electronics and the antennas at a propagation speed near that of the wave speed 
in the host medium. Accordingly, unless compensated for, this can prevent accurate absolute time 
registration, i.e., determining the actual time at which the transmitted wave leaves the transmitting 
antenna. Moreover, in an array-based system, electronic circuitry can differ for each antenna pairing. 
Consequently, an embodiment of data processor 1 04 consistent with the present invention performs 
time alignment for each antenna pairing used in sensor unit 102. Moreover an embodiment of data 
processor 104 consistent with the present invention can include a wave speed calculator to compute 
wave speed as described in greater detail below. 

11 
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[039] Time alignment is accomplished by recognizing that the wave sent from each 
transmitting antenna to each receiving antenna travels through air. The distance between any 
transmitter-receiver antenna pair is know. In Fig. 8 , the distance between transmitting antenna 
element 802 and receiving antenna element 804 is labeled "d. " The speed of light in air is also 
know: c air = 3 x 10 8 meters per second. Accordingly, data processor 104 computes the time of this 
air propagated arrival using the following equation: 



(1) 



Because the wave speed in air is always greater than or equal to the speed of light in geological 
material, the earliest arriving signal will be from the air-propagated wave. For each transmitter- 
receiver pair, data processor 1 04 searches an acquired time series (Fig. 9) for the first arriving peak. 
As seen in Fig. 9, the peak at approximately twelve nanoseconds is the air-propagated arrival time 
and the peak at approximately forty nanoseconds in the ground-propagated arrival time. Data 
processor 104 then adjusts the acquired time series so that the arrival time of the first peak occurs 
at t air as determined from Eq. (1). 

[040] A GPR apparatus consistent with the present invention uses spatial correlations 
among time series acquired for various transmitter-receiver antenna element pairings. The travel 
distance r is related to the travel time t arrival by: 



r ~ C <f arrival > 



(2) 
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where c 0 is the wave speed in the host medium. Data processor 1 04 computes the background wave 
speed c 0 . 

[041] Data processor 104 computes c 0 by exploiting the fact that there will be a critically 
refracted ray (Fig. 1 0, 1 008) that will propagate at speed c„ horizontally between a transmitter ( 1 002) 
and receiver (1004) just below the air-ground interface (1006). This ground-propagated arrival will 
occur later than the direct air-propagated arrival. Since this is the shortest ground-propagated ray 
path, this arrival will precede any ray paths associated with reflections from subsurface objects. This 
concept is illustrated in Fig. 11. Fig. 11 shows critically refracted ray path 1108 and ray path 
associated with a reflection 1 1 10a, 1 1 10b from subsurface object 1 106. Since ray path 1 108 for the 
critically refracted ray is shorter than the ray path associated with a reflection 1 1 10a, 1 1 10b, the 
arrival of the ground-propagated ray will precede that of the reflection. 

[042] Fig. 1 2 is a schematic diagram of bistatic measurements for computing wave speed. 
The distance from the transmitter to receiver i is denoted by d,. For example, in Fig. 12, d, is the 
distance 1206 between transmitting antenna (T) 1202 and receiving antenna (R, ) 1204; d 2 is the 
distance 1210 between transmitting antenna (T) 1202 and receiving antenna (R 2 ) 1208; and </„isthe 
distance 1214 between transmitting antenna (T) 1202 and receiving antenna (R„ ) 1212. Note also 
that when a transmit/receive antenna pair is used instead of an array as described above, the distance 
used is the distance from the transmitter to the receiver position. 

[043] Data processor 1 04 computes the arrival time for the critically refracted ray at receiver 
i using the equation: 




(3) 
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As previously noted, the first major arrival is the air propagated arrival and the second major arrival 
is the critically refracted ground propagated arrival. Data processor 104 searches the time series for 
all receivers to find the t,'s. If time alignment has been performed, data processor 104 can use Eq. 
(3) to compute c 0 directly by selecting a particular receiving antenna and, with /,. and d t known using 
the following equation: 



(4) 



t. 



Alternatively, all f/s and d/s can be used and — is the slope of a linear fit for Eq. (3), e.g., the best 



linear fit. 



[044] If data processor 1 04 has not performed time alignment, a similar procedure can be 
applied using the time difference At between the air and ground propagated arrivals. In this case, 
the appropriate relationship is given by the equation: 



a d i d i 
A/. = — - — 



C 0 C air 



(5) 



rather than Eq. (3). 

[045] When using an array of antennas, there may be arrivals intermediate to the air- and 
ground-propagated arrivals. These are typically air-propagated and associated with reflections or 
multiple reflections with the antenna array. These artifacts are referred to herein as "ringing" and 
will be manifested as peaks in the power spectra of the acquired data. It may be necessary to remove 
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this ringing before computing wave speed. This can be done by identifying the spectral peaks and 
applying appropriate spectral band-pass filters, such as Blackman-Harris filters, known to those 
skilled in the art prior to implementing any of the above-described procedures. 

[046] The following is an additional method for computing wave speed consistent with the 
present invention in a GPR system configured with a single transmitting antenna and one or more 
receiving antennas where the i th receiver is a distance d s away from the transmitting antenna (or when 
one receive element is used and moved to simulate an array, the distance used is the distance from 
the transmitter to the receiver position). If the transmitter is pulsed, the first two arriving signals at 
any receiver are transmitted through the air propagating at the wave speed in air, c Q , and a critically 
refracted ray propagating just below the air-ground interface at a wave speed c g , where c g *c a . 

[047] Knowing the distance between the transmitter and receiver, the wave speed in the 
ground can be computed using the equation: 



g z n 



(6) 



where d is the transmitter-receiver separation and x g is the measured travel time of the critically 
refracted ray through the ground. This procedure is useful when, for example, 1) the peak in 
amplitude as a function of time associated with the critically refracted ray can be clearly identified 
and 2) there is time registry in the acquired data. For radio waves, the information is traveling 
through the GPR system circuitry at about the same speed as the wave, time registry can be difficult 
resulting in an error in x g that will introduce an error in c g (Eq. 6). Since a time registry problem 
produces a time shift to the entire time series, the travel time x g will be effected; however, the 
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difference in travel time between the critically refracted ray and the air-propagated ray t g . fl = x g - x a 
will not be effected by this error. Consequently, data processor 104 uses this time difference to 
compute the wave speed c g using the equation: 



(7) 



where c a and d are know so that 



c = 



X + I 

g-a a 



(8) 



[048] When there is substantial noise in the data (e.g., condition 2, above), it may not be 
possible to definitively identify the appropriate peaks within the time series so that direct 
implementation of Eq. (8) is impossible. In such cases, can be determined indirectly. Let Stf) 
be a time series for a transmitter-receiver spacing of d g and define s^co) to be the temporal Fourier 
transform of S { , 



(9) 



An auto-correlation of S ; can have the form 



e 



-KD/ 



j,(<oK(co) 
j,.((o K(to)+e 



(10) 



where the asterisk denotes complex conjugation and e is small parameter to eliminate singilarities 
and supress noise in the data. A plot of the time series F t can have three peaks; one 
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at t = 0 corresponding to the correlation of the air propagated arrival with itself and the ground 
propagated arrival with itself; as well as peaks at ± (x a - x g ) corresponding to the correlations 
between air and ground propagated peaks. These peaks can be used to determine x g ^ and Eq. (8) can 
then be employed to compute c g . 

[049] An embodiment of data processor 104 can implement a variation of the above 
procedure using two receiving antennas at different distances from the transmitting antenna. Let 5, 
and S 2 be time series for transmitter-receiver spacings of d x and d 2 , respectively. A cross-correlation 
of the form 



J, (©)$,*(©)+€ ' 



(11) 



will have significant peaks at \\ - x 2 a ,x\ - - %\ , and t^-t^, where T' a and x' g are 

the air and ground propagated travel times, respectively, for the i' h receiver. These time differences 
can be used by data processor 1 04 in various ways to compute c g . For example, using the difference 

t g - t g ; because d x = c g x g and, data processor 104 computes the wave speed using the following 
equation: 



d x -d 2 

8 



(12) 



[050] Another calibration operation that can be performed by data processor 104 is the 
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computation of an impulse response function for the GPR system. Ground penetrating radar 
antennas have spatially three-dimensional beam patterns that are frequency dependent. While it is 
possible to characterize these beam patterns by direct propagation measurements in air, these patterns 
will change when the antennas are dynamically coupled to a material with a dielectric or electrical 
conductivity that is different from air. In an embodiment of the present invention, data processor 1 04 
characterizes antenna beam patterns from any array geometry by acquiring data over an obj ect having 
a known location and known electromagnetic properties. 

[05 1 ] In the frequency domain, the acquired data at some point r in the array of receivers 
for a transmitter position r„ is expressed as u (r, r 0 ; where k 0 = co/c 0 is the background 
wavenumber at frequency <o. The wave field u can be considered as the superposition of u 0 , the 
wavefield that would exist in the absence of a buried object and the perturbation to u 0 associated 
with the known buried object. With these definitions, the perturbed wave field is: 



«'(r, r 0 ; *„) = - ^ J dr' G(r - r'; A>(r'-r 0 ; k 0 )O(r), 



(13) 



where G is the Green's function and O is the "object function" that characterizes the buried object: 



0(r) = 1 



(14) 



* 2 00 = M 0 e(r)co 2 + ifi<p(r)a), 



(15) 



In this equation: ^ 0 is the free-space magnetic permeability, e 0 and e(r) are background and object 
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dielectric constants, respectively; and a Q and o(r) are the background and object electrical 
conductivities, respectively. Since any object can be buried, a relatively small object can be selected 
such that the Born approximation is valid, in which case, the background wave field u 0 can be used 
in place of the total field u in Eq. (1 3) to give 



u'(r, r 0 ; k 0 ) = - ^ J dr 1 G(r-r'; *> 0 (r ' - r 0 ; k 0 )O(r). 



(16) 



[052] Equation (16) can be used by data processor 104 for forwarding modeling where, 
given 0, u' is computed. Alternatively, data processor 104 can use this equation for inverse 
scattering where O is determined from known u 0 and measured u\ Here, it is assumed that u 0 and 
O are known, u' is measured for a measurement configuration of interest, and Eq. (16) is used by 
data processor 104 to determine the impulse response function of the antenna array that includes the 
effects of the unknown beam pattern. 

[053] For point (monopole) receivers, the Green's function can be represented by the Weyl 
expansion 



(17) 



where z is a unit vector in the vertical direction. For any arbitrary receiving antenna beam pattern, 
the Green's function can be written as 



G(r-r'; k 0 ) = i j — — — A R (a; k y^^ 1 '^ 



where the factor A R is included to represent the unknown receiving antenna beam pattern. 
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[054] Similarly, the transmitting antenna beam pattern can be introduced into the 
background field u 0 in the form 



Wo (r'-r 0 ; ^ = ij -^—A T (fi; Q.V+'^-fiTW-* 



(19) 



Substitution of Eqs. (18) and (19) into Eq. (16) yields a complicated relationship between u' and 
O involving the product A^Aj, For a fully bistatic geometry, the horizontal two-dimensional Fourier 
transform 



ff'Or. v; * 0 ) = j dr j dr 0 M '(r, r 0 ; k^'™ •> 



will "deconvolve" Eq. (16) yielding the linear relationship 



u'(k, v; k 0 ) = A(k, v; k Q )I(( K , v; k 0 )O, 



(20) 



(21) 



where i"is the two-dimensional spatial Fourier transform of AgA^ lis the impulse response function 
for monopole sources, and 6 is the three-dimensional spatial Fourier transform of the object 
function. 

[055] Accordingly, data processor 104 uses Eq. (21) to compute the spatial Fourier 
transform of the impulse response function: 



u'{k, v; k Q ) 



(22) 
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for the array including the antenna beam patterns and spectral respnses in terms of the known O and 
measured i/'. 

[056] Arrays that are not fully bistatic are either fully monostatic or bistatic in one direction 
and monostatic in the other. Accordingly, for these geometries, additional assumptions or 
approximations can be used by data processor 104 for deconvolving Eq. (16). One suitable 
assumption is that the object is buried about one wavelength deep or deeper. Since the burial depth 
of the reference object can be controlled during calibration, the validity of this assumption can be 
assured and a valid relationship such as that given by Eq. (16) can be developed. Thus, using the 
analytical framework described above, data processor 1 04 detects spatial correlations corresponding 
to the buried object in data from the GPR sensors. 

[057] The next element of the GPR systems of Fig. 1 is preprocessor 106. An embodiment 
of preprocessor 106 consistent with the present invention is comprised of a processor that executes 
stored computer program code designed to implement preprocessing operations. One skilled in the 
art will recognize that an embodiment of preprocessor 1 06 implemented entirely in software, entirely 
in hardware, or in an embodiment allocating preprocessing functions among hardware and software 
elements, either distributed or centralized, is consistent with the scope of the present invention. 

[058] Preprocessor 106 rejects elements in the acquired data associated with unwanted 
subsurface features and/or admits only targets of interest. For example, consider spatial frequency 
preprocessing of monostatic data. Frequently reflections from soil or geological strata produce 
strong reflections that can obscure other features of interest. These strata typically are nearly planar 
and horizontal. The two-dimensional spatial Fourier transform of a horizontal surface will be a delta 
function at the spatial frequency origin. An embodiment of preprocessor 106 consistent with the 
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present invention removes such a feature using a notch filter for filtering a DC response. Since most 
naturally occurring strata are not perfectly horizontal, and may not be truly planar, preprocessor 1 06 
relaxes the notched filter somewhat into a low-cut filter (e.g., "high-pass filter) to substantially 
suppress these strata. A Blackman-Harris filter is an example of an appropriate filter for this 
purpose. Similar spatial frequency filters can be used to admit features of interest. If targets of 
interest are buried utilities, for example, these have the rather unique characteristic of being long and 
thin. The spatial Fourier transform of such a shape will generate a DC response in the axial direction 
and a white response in the transverse direction. Thus, preprocessor 1 06 applies an appropriate filter 
to admit features having these characteristics. 

[059] Image processor 1 08 builds a data structure corresponding to an image of the buried 
object from the data processed by data processor 104. Image processor 108 can, for example, 
reconstruct the location, size, shape, and material properties of subsurface features in two- 
dimensions, three-dimensions, or as a sequence of two-dimensional images. An embodiment of 
image processor 1 08 consistent with the present invention is comprised of a processor that executes 
stored computer program code designed to implement image processing operations. One skilled in 
the art will recognize that an embodiment of image processor 1 08 implemented entirely in software, 
entirely in hardware, or in an embodiment allocating image processing functions among hardware 
and software elements, either distributed or centralized, is consistent with the scope of the present 
invention. 

[060] Embodiments of image processor 108 implement image processing algorithms 
including, but not limited to, synthetic aperture radar (SAR), migration, backpropagation, diffraction 
tomography, algebraic reconstruction tomography (ART), iterative ART, pseudo-inverse, and reverse 
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time processing. One skilled in the art will recognize that other image processing algorithms are also 
suitable for implementing the functions performed by image processor 108. Moreover, some 
operations performed by preprocessor 106 can be incorporated directly into many of these imaging 
procedures of image processor 108 by introducing filters into the impulse response function. 

[061] Feature processor 110 processes a received GPR signal to characterize features 
present in the signal that correspond to properties of a detected buried object. An embodiment of 
feature processor 1 10 consistent with the present invention is comprised of a processor that executes 
stored computer program code designed to implement feature processor operations. One skilled in 
the art will recognize that an embodiment of feature processor 1 1 0 implemented entirely in software, 
entirely in hardware, or in an embodiment allocating feature processor functions among hardware 
and software elements, either distributed or centralized, is consistent with the scope of the present 
invention. 

[062] Feature processor 1 10 differentiates the GPR signal to characterize features based 
on differences in reconstructed material properties. An embodiment of feature processor 1 10 uses 
frequency domain processing, performing an inverse operation mapping frequency information to 
depth. When feature processor 110 processes a GPR signal that propagated through a non- 
conducting background, the reconstructed quantity is the object function, Eq. (14), which has the 
form 



0( X) = i +< £&. 



(23) 



e 0 coe 0 



In Eq. (23), O is complex and its real part, 
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■ e(x) 



(24) 



is a reconstruction of the spatial variations in dielectric constant. The imaginary part of 0, 
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explicitly involves the frequency. 

[063] Let L" 1 be the operator that maps the broadband data into the image, i.e., 

0 = L l {u'} 

is the inverse operator, so that L is the forward operator, 



k' = L{0}. 



Using Eq. (23): 



where 



O a (x) = V l {coLicoQ,}}, 
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and the pair of images O e ,O a characterize imaged features by their relative spatial variations in 
dielectric and electrical conductivity. Those skilled in the art will recognize that equations 23-27 
can be extended to accommodate processing a GPR signal that propagated through a conducting 
background. 

[064] Feature processor 1 1 0 accomplishes further differentiation by recognizing how O e and 
O a are related for certain material types. For example, plastic objects are characterized by low 
dielectric and low conductivity. Consequently, feature processor 1 1 0 recognizes that plastic objects 
exhibit an O e that is greater than background (positive) and an O a that is, in general, less than 
background (negative). The converse is generally true for metal objects. For this reason, feature 
processor 110 generates an enhanced data set forming a picture of metal or plastic objects by 
examining O e - For metal objects, this difference could be strongly negative while, for plastic 
objects, it could be strongly positive. For other materials such as clay, wood, and concrete; feature 
processor 110 enhances features using the sum O e + O a . 

[065] Detector 1 1 2 automatically identifies targets of interest. An embodiment of detector 
1 12 consistent with the present invention is comprised of a processor that executes stored computer 
program code designed to implement detection operations. One skilled in the art will recognize that 
an embodiment of detector 1 12 implemented entirely in software, entirely in hardware, or in an 
embodiment allocating detection functions among hardware and software elements, either distributed 
or centralized, is consistent with the scope of the present invention. 

[066] Detector 1 12 determines whether an object detected in a surveyed region matches a 
predetermined target or belongs to a predefined target class. If the object in the surveyed region does 
match, detector 112 identifies the location of the object. An embodiment of detector 112 uses 

25 



Fl NNEGAN 
HENDERSON 
FAR A BOW 
CARRE TT & 
DUNNER lip 

13001 Street, NW 
Washington, DC 20005 
202.408.4000 
Fax 202.408.4400 
www.finnegan.com 



o 



o 



matched filters to compare a signal representing an object in a surveyed region to targets and target 
classes. In one embodiment of detector 1 12 using matched filters w(r, co) is the acquired data over 
a measurement surface defined by r over a bandwidth spanned by co. Similarly, u(r, co; x*) is the 
data that would be acquired if the target of interest were centered at the point Xq. This target is 
present in the measured data set if 



F Mo) = £ J dru (r, ©Mr, co; x 0 ) 



(28) 



exceeds a threshold and the location of this target will be at the position Xq that maximizes F c . 

[067] Many other detection algorithms are suitable for implementation by detector 1 12 
such, for example, Baysian statistics or the MUSIC algorithm. 

[068] Descriminator 1 14 is a signal processing element similar to detector 1 12 (described 
in detail above); however, descriminator 114 processes images reconstructed from received GPR 
signal data rather than processing the acquired GPR signal data directly. An embodiment of 
descriminator 1 14 consistent with the present invention is comprised of a processor that executes 
stored computer program code designed to implement descriminator operations. One skilled in the 
art will recognize that an embodiment of descriminator 114 implemented entirely in software, 
entirely in hardware, or in an embodiment allocating descriminator functions among hardware and 
software elements, either distributed or centralized, is consistent with the scope of the present 
invention. 

[069] There are a number of approaches to discrimination consistent with the present 
invention. For example, an embodiment of descriminator 114 performs a correlation operation 
similar to Eq. (28), but applied to images reconstructed from received GPR signal data. Other 
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embodiments of descriminator 114 implement feature extraction algorithms, pattern recognition 
algorithms, and/or feature tracking algorithms that track a feature identified in one portion of the 
image throughout the entire image. 

[070] Display 1 1 6 proves an operator of the GPR system with visual feedback representing 
the surveyed region. An embodiment of display 116 consistent with the present invention is 
comprised of a graphical display device that can present an image to an operator (e.g., video screen, 
plotter, printer, etc.) and a processor that executes stored computer program code designed to 
implement display operations. One skilled in the art will recognize that an embodiment of display 
1 16 implemented entirely in software, entirely in hardware, or in an embodiment allocating display 
functions among hardware and software elements, either distributed or centralized, is consistent with 
the scope of the present invention. 

[07 1 ] In an embodiment of display 1 1 6, the image presented to the operator takes the form 
of a three-dimensional rendering of the surveyed site or a sequence of two-dimensional images of 
the reconstruction. The reconstructions can be of a basic material property, such as dielectric, or a 
composite of multiple properties, such as reflectivity, object function or some statistical estimator. 
Other embodiments of display 116 consistent with the present invention present images to the 
operator as projections of reconstructed information onto a plane such as the ground surface. 

[072] The GPR system functions described for each of the elements above can be 
distributed or combined differently than stated above and still fall within the scope of the invention. 
For example, the GPR system can be implemented using a single processor or group of processors 
that execute software stored in a central memory, wherein the software relates to each of the above- 
described GPR system operations. 
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[073] The operation of the GPR system elements as defined above provide an efficient 
system to rapidly and accurately locate buried objects. Therefore, in accordance with the present 
invention, the disclosed GPR system elements can be housed in many types of structures, including 
portable and mobile structures to allow rapid and accurate real-time surveying of a region of interest. 

[074] Fig. 1 3 is a diagram of a mobile vehicle 1 304 with a GPR system, consistent with this 
invention. A radar array 1306 attaches to an arm 1305, which attaches to the back a vehicle 1304, 
as shown in Fig. 13. Radar array 1306 may comprise at least one transmit antenna and at least one 
receive antenna. Vehicle 1304 may move in direction shown by arrow 1302. Radar array 1306 
transmits impulses into ground 1308. The impulses may reflect off of a buried object, such as 
subterranean pipe 1 3 1 2 and radar array 1 306 may receive reflected waveforms. Module 1 3 14 on the 
back of vehicle 1304 may comprise elements described in greater detail above for Fig. 1, which 
control array 1306 and process signals received by array 1306. It may also display on display 1316 
images of buried objects from data acquired by sensor unit 102 in the surveyed region of interest 
(with or without images of the surrounding area). 

[075] Fig. 14 is a diagram, consistent with this invention, of mobile vehicle 1404 with a 
trailer 1 402 having a GPR system. In the embodiment shown in Fig. 1 4, radar array 1 306 (not shown 
in Fig. 14) is within trailer 1402, which is attached to vehicle 1404. Trailer 1402 moves in the 
direction of arrow 1410 with vehicle 1404. Fig. 15 is a diagram, consistent with this invention, of 
a portable housing 1504 with a GPR system. In the embodiment shown in Fig. 1 5, radar array 1306 
(not shown in Fig. 14) is within portable housing 1504. A user 1502 may guide portable housing 
1504 over ground 1308 using handle 1506. Portable housing 1504 may have wheels 1508, skid pads 
(not shown), or neither. An embodiment of portable housing 1 504 can also be made sufficiently light 
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to allow an operator to carry the entire system, thereby alleviating the need for wheels, skid pads, 
similar structure to guide the housing over the ground. Both the mobile system of Fig. 14 and the 
portable system of Fig. 15 may also display images of buried objects from data acquired by the 
sensor unit in the surveyed region of interest (with or without images of the surrounding area). 

[076] A GPR system consistent with the present invention can overcome the disadvantages 
of One Call, private locating, and SUE by providing a cost effective method to locate and image 
conductive and non-conductive utilities, vertically and horizontally, with a margin of error to satisfy 
many excavating needs. A GPR system consistent with the present invention can also provide a 
permanent record of images of the excavation site that can be used in the future. Further, a GPR 
system consistent with the present invention can fit onto a moving vehicle, trailer, or portable 
housing so that subsurface images can be formed as the system is moving. 

[077] It will be apparent to those skilled in the art that various modifications and variations 
can be made to the disclosed embodiments of the present invention without departing from the spirit 
or scope of the invention. Thus it is intended that the present invention cover the modifications and 
variations of this invention provided they come within the scope of the appended claims and their 
equivalents. 
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